Interacting anyons in topological quantum liquids: The golden chain 
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We discuss generalizations of quantum spin Hamiltonians using anyonic degrees of freedom. Tiie simplest 
model for interacting anyons energetically favors neighboring anyons to fuse into the trivial ('identity') channel, 
similar to the quantum Heisenberg model favoring neighboring spins to form spin singlets. Numerical simula- 
tions of a chain of Fibonacci anyons show that the model is critical with a dynamical critical exponent z = 1, 
and described by a two-dimensional (2D) conformal field theory with central charge c — 7/10. An exact map- 
ping of the anyonic chain onto the 2D tricritical Ising model is given using the restricted-solid-on-solid (RSOS) 
representation of the Temperley-Lieb algebra. The gaplessness of the chain is shown to have topological origin. 

PACS numbers: 05.30.Pr, 73.43.Lp, 03.65.Vf 



Introduction Non-Abelian anyons are exotic particles ex- 
pected to exist in certain fractional quantum Hall (FQH) states 
2J. A set of several anyons supports very robust collec- 
tive states that are degenerate to exponential precision; such 
states can potentially be used as quantum memory and for 
quantum computation However, this degeneracy can be 
lifted by a short-range interaction if the anyons are very close 
to each other. As a first step towards understanding interacting 
anyons, we describe a simple, exactly solvable model that is 
an anyonic analogue of the quantum Heisenberg chain. 

We start by considering the well-known Moore-Read state 
HI], a candidate state, exhibiting non-Abelian statistics, for the 
topological nature of FQH liquids at filling fraction u = 5/2. 
It has two important types of excitations: quasiholes with 
electric charge e/4 and neutral fermions. Quasiholes may be 
trapped by an impurity potential while the fermions can still 
tunnel between them |4]. For a one-dimensional (ID) array of 
trapped quasiholes, the Hamiltonian can be described in terms 
of free Majorana fermions on a lattice, which is in turn equiv- 
alent to the ID transverse field Ising model at the quantum 
phase transition point. The more interesting model discussed 
here is based on so-called 'Fibonacci anyons', which repre- 
sent the non-Abelian part of the quasiparticle statistics in the 
fc = 3, Zfe-parafermion state [2^, an effective theory for FQH 
liquids at filling fraction v = 12/5 [jst]. Even without pa- 
rameter fine-tuning, these ID anyonic arrays will be shown to 
exhibit gapless excitations due to topological symmetry. 

Model Our model describes pairwise interactions within 
an aiTay of L anyons, for instance along a chain as shown in 
Fig. [T^). In the Fibonacci theory there are only two types of 
particles: the Fibonacci anyon, denoted by r, and the trivial 
particle denoted by 1 with a fusion rule r x r = 1 + r. We 
refer to the label 1 or r as the topological charge. When two 
neighboring anyons interact, indicated in the figure by the el- 
lipses, they can either fuse in the trivial channel, annihilating 
each other, or in the nontrivial one, becoming a single r-anyon 



f^. We define our model by assigning an energy gain if they 
fuse along the trivial channel. This is an anyonic analogue of 
the spin- 1/2 quantum Heisenberg antiferromagnet, which as- 
signs an energy gain to two neighboring spin- 1/2 fusing into a 
spin-0 singlet as compared to a spin- 1 triplet. 

To define the Hilbert space of r-anyons we consider the 
tree-like fusion diagram in Fig. [TJ)). The basis corresponds to 
all admissible labelings \xi, X2, . . .) of the links, with Xi — 1 
or r. Each label represents the combined topological charge 
of the particles left to a given point. Not all possible values 
{xi,X2, . . .) represent allowed basis states due to the fusion 
rules: a 1 must always be preceded and followed by a r, since 
the fusion of a 1 and a r always gives a t. This reduces the di- 
mension of the Hilbert space of the open chain (with r-labels 
at the boundary) to the Fibonnacci sequence dim^ = -Fl+i, 
and for the periodic chain dim^ = Fl-i + Fl+i- For 
large L it is well-known that these numbers grow at a rate 
dim^ oc i/j-^, where ip = {1 + V5)/2 is the golden ratio. 
This Hilbert space has no natural decomposition in the form 
of a tensor product of single-site states, in contrast to SU (2) 
quantum spin chains. 
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FIG. 1: a) Illustration of the Fibonacci chain with L r-anyons. 
b) The fusion path, c) Definition of the F-matrix. 
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In order to generate a local Hamiltonian H = ^ ■ H; as- 
signing an energy to the fusion of two neighboring r-anyons 
we use the so-called F-matrix to transform the local basis as 
shown in Fig. [T];). In the transformed basis the state x'^ cor- 
responds to the fusion of the two anyons. The Hamiltonian is 
then defined by assigning an energy Er — for x'^ = t, and 
El = —1 for x[ = 1. The resulting local terms H,; contain 
three-body interactions in the link basis, 

Hilxi-ix^Xi+i) = ^ (H.i)x," \xi-ix' iXi+i) (1) 

X' i—l.T 

with (H^)-; HF^:^:rr)liF^:^:rrt^ ■ 

It is diagonal in the subspace {\xi-iXiXi+i)} = 
{\1t1),\1tt),\tt1)}. Hi = diag{-l,0,0}, where the 
i^-matrix is a number due to the constraints arising from 
the fusion rules. For the case xi-i = Xj+i = r, the F- 
matrix and the corresponding Hamiltonian are the following 
2 X 2-matrices (x^, e {1,t}) 

(2) 

Looking at the matrix form of the Hamiltonian, it can be writ- 
ten in terms of standard Pauli matrices: 

H, = + - 1) 

-ni^im+i {ip^'^l'^al + ip-'^rii + 1 + (^"^^ , 

where the sum runs over the links of the chain. In this expres- 
sion, the operators Ui count the r-particle occupation on link 
i, Tij = i(l — erf) = 0, 1, and the Hamiltonian H acts on the 
constrained Hilbert space defined above. 

Entropy scaling and central charge We simulated this 
model numerically, and calculated the finite-size excitation 
gap by exact diagonalization for chains of up to L = 32 sites. 
A finite-size analysis shows that the gap vanishes linearly in 
1 /L, indicative of a critical model with dynamical critical ex- 
ponent z = \ described by a conformal field theory. The 
central charge c of a CFT can be calculated from the finite- 
size scaling of the entanglement entropy. For two subsystems 
with equal size L/2 on systems with periodic (PBC) and open 
boundary conditions (OBC) the entanglement entropy scales 
as 101 

c c 
5'pBc(i) oc - logL and S'oBc(i) oc - logL . (3) 
3 6 

The density matrix renormalization group method (DMRG) 
lUlSl provides a natural framework for calculating these quan- 
tities. Fits of our numerical results according to Eqs. (|3]l 
shown in Fig. |2] give central charge estimates of cpBc = 
0.701 ± 0.001 and cqbc = 0.70 ± 0.01 respectively. Since 
possible (unitary) CFTs in the vicinity of these estimates 
have central charges|19] 1/2, 7/10 or 4/5 we can unambigu- 
ously conclude that our results are consistent only with central 
charge c = 7/10. 



Mapping and exact solution We now proceed to derive 
these results exactly. By construction the local contribution to 
the Hamiltonian ^Xj = — is a projector onto the trivial 
particle. One can then verify that the operators form a 
representation of the Temperley-Lieb algebra I1IQ1I 

(xo' = dx„ 

X^XiiiXj ~ Xj, 

[X„X,] = for > 2, (4) 

where the 'd-isotopy' parameter equals the golden ratio, d — 
(p. This representation can be seen to be identical to the stan- 
dard Temperley-Lieb algebra representation associated with 
SU{2)k at level fc = 3. For arbitrary k, the latter contains 
k + 1 anyon species labelled by j = 0, 1/2, 1, A:/2, satisfy- 
ing the fusion rules of SU{2)k 1 1 1]. The operators defined 
by 



and 




are known lfl2ll to form a representation of the Temperley-Lieb 
algebra for any value of k, where \ji — = 1/2 and 

Our model of interacting Fibonacci anyons can be cast into 
this form at fc = 3 by first mapping xi = 1 ^ ji = 0, and 
Xi — T ji — 1, and then applying the SU{2)j, fusion rule 
3/2 X j = 3/2 — j to the even-numbered sites. This maps any 
admissible labeling \x) := \xi,X2, ■ ■ ■) uniquely into |j) := 
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FIG. 2: (Color online) Entropy scaling for interacting Fibonacci 
anyons arranged along an open (open squares) or periodic chain 
(closed circles) versus the system size L. Logarithmic fits (solid 
lines) give central charge estimates of cpBC = 0.701 ± 0.001 and 
coBC = 0.70 ± 0.01 respectively, where for the open boundary con- 
ditions only the values for the 5 largest systems have been taken into 
account due to large finite-size effects. 
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I Ji , ^2 , • ■ ■) where for odd-numbered sites j2i+i G {0,1}, and 
for even numbered-sites G {1/2,3/2}. This re-labeling 
maps the matrix elements of Xj into those of from Eq. (|5]). 

We can now see that the Hamiltonian in Eq. ([T]i is that cor- 
responding to a standard (integrable) lattice model descrip- 
tion of the classical 2D tricritical Ising model, known as the 
RSOS model t^- Specifically, the two-row transfer matrix 
T := T2T1 of this lattice model, shown in Fig. [3] is written 
in terms of Boltzmann weights W[i] assigned to a plaquette i 
of the square lattice 

Ti:=]JW[2n], and T2 := J| W[2n + 1] 



with 



T sm T-frr — u -, sin?/ 

W[i]i = \t + eMi 

' sin^ J sin^ J- 



(6) 



The parameter m > is a measure of the lattice anisotropy, 
1 is the identity operator, and 



The Hamiltonian of the so-defined lattice model is obtained 
from its transfer matrix by taking, as usual l,15i] . the extremely 
anisotropic limit, u <C 1, 



T = exp{-a(H + ci) + 0(0^)}, a = 



uip 



sin[7r/(fc + 2)] 



< 1 



yielding H = — ^ - ^e^ (ci is an unimportant constant). 
Since the operators X; can be identified with e;, this demon- 
strates that the Hamiltonian of the Fibonacci chain is exactly 
that of the corresponding k = 3 RSOS model which is a lattice 
description of the tricritical Ising model at its critical point. 
The latter is a well-known (supersymmetric) CFT with cen- 
tral charge c ~ 7/10 131 ■ Analogously one obtains 1 181 
for general k the (fc — 1)^* unitary minimal CFT 1 19] of central 
chargec ~ 1 — 6/(fc + l)(A; + 2). A fermmagnetically coupled 
Fibonacci chain (energetically favoring the fusion along the r- 
channel) is described by the critical 3-state Potts model with 
c = 4/5 and, for general k, by the critical Z^-parafermion 
CFT lUlIillli with central charge c = 2(fc - l)/{k + 2). 




FIG. 3: Transfer matrix of the RSOS model. 
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FIG. 4: Energy spectra for periodic Fibonacci chains of size L = 36 
and L — 37. The spectra have been rescaled and shifted such that 
the two lowest eigenvalues match the conformal field theory assign- 
ments. The open boxes indicate the positions of the primary fields 
of the c — 7/10 conformal field theory. The open circles give the 
positions of multiple descendant fields as indicated. While we find 
excellent agreement in general, finite-size effects lead to small dis- 
crepancies for the higher energy states. The solid line is a cosine-fit 
of the dispersion which serves as a guide to the eye. 



Excitation spectra We have calculated the excitation 
spectra of chains up to size L = 37 with open and periodic 
boundary conditions using exact diagonalization, as shown in 
Fig. E] The numerical results not only confirm the CFT pre- 
dictions but also reveal some important details about the cor- 
respondence between continuous fields and microscopic ob- 
servables. In general, low-energy states on a ring are associ- 
ated with local conformal fields |21], whose holomorphic and 
antiholomorphic parts belong to representations of the Vira- 
soro algebra, described by conformal weights and hj^. The 
energy levels are given by 
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symmetries. Such terms are relevant if h'f^ + < 2, in 
which case they may open a spectral gap or induce crossover 
to different critical behavior at large distances. In the tricriti- 
cal Ising model, there are four relevant fields: £l£_r, e' l^' r, 
ctlO'r, l<^' r- Some explanation is in order as to why these 
fields do not appear in the effective Lagrangian of our model. 
The fields a and a' are staggered and thus prohibited by trans- 
lational symmetry. Excluding e and e' requires a more subtle 
argument. The Fibonacci ring has a topological symmetry, 
which corresponds to adding an extra r-line parallel to the 
spine of the fusion diagram (Fig.[TJ)) and merging it with the 
diagram using the i^-matrix. We denote this operator by Y . 



TABLE I: Lowest eigenvalues for open Fibonacci chains of size L = 
31 and L = 32. The two lowest eigenvalues have been rescaled and 
shifted such that they match the conformal field theory assignments. 



corresponding to states with a choice of momenta K — Hl - 
hji or K = hi — hj^ + L/2 in units of 2tt/L, where Ei 
V are non-universal constants. Here, — 



h^?^ + rriL and 



where h^j^^ , /i^"* correspond to weights of 
'primary' fields and and are non-negative integers 
describing so-called 'descendant' fields. The numerical spec- 
tra for even values of L (see the first plot in Fig. IH agree 



with Eq. (I8]l, exhibiting primary fields with h)^ = h]^ = 
0,1/10,3/5,3/2,3/80,7/16, which are conventionally de- 
noted by /, e, e', e", cr, a', respectively [24J. The momenta of 
the last two fields and their descendants are near K = L/2, 
as compared to the other four, indicating that the correspond- 
ing microscopic observables have alternating sign on the lat- 
tice. Such "staggered" fields must have nontrivial monodromy 
with respect to a space-time dislocation (i.e., the insertion or 
removal of a site at some particular time). Such a dislocation 
is characterized by a chiral e" field, say, e'l [22]. The log of 

the monodromy factor exp(27ri (/i^ ^ ^ ~ ^11" )) matches 
the momenta K in Fig. |1] f23\. Given this information, we 
may predict that the states of an odd size ring are associated 
with fields of the form ^lVR' where 77 = ^ x e". These include 
six primary fields, e'/, e^, e'j^Ej^, Sj^e'j^, cr^o-j^, cr^cr^, as well 
as their descendants. Integrality of the momentum K dictates 
the choice of K (see below Eq. (HJ), as in Fig.|4] 

For open boundary conditions the spectra are known to be 
described by, say, the holomorphic sector only [26J. To ex- 
plain the numerical data, we need to assume that the ends of 
the chain are charaterized by a boundary field e (or equally 
well e'). Thus, for an even number of sites the spectrum is 
described by e x e = / + e' (plus descendants). For an odd 
number of sites, this result is to be modified by fusion with 
e", yielding e" + e. The numerical low-energy spectra agree 
excellently with these predictions as shown in Table |T] 

Hidden symmetries The critical behavior of our model is 
not just a peculiarity of the exact solution but rather has topo- 
logical origin. In general, an effective low-energy Lagrangian 
admits perturbations of the form J ipix^r) dx dr, where -if) 
may be any local field that is consistent with all applicable 



L-l 



1=0 



where the identification L = is used. We may think of the 
fusion diagram as a description of a process that generates a 
set of r-anyons on a circle from the local vacuum. Then Y 
describes another particle moving inside or outside that circle 
(or on the circle itself — before those anyons were created). 
The operator Y is sensitive to a possible topological charge 
y = 1,T located at the center of the circle. Thus Y has two 
eigenvalues, Syr/Syi = ip,—ifi^^. We conjecture that the 
low-energy states associated with fields /, e", cr' are in the the 
trivial (y — 1) sector, and the fields e',e,a are in the y = t 
sector. In fact, the topological fusion algebra (defined by the 
rule TXr = l + T)isa quotient of the CFT fusion algebra. 

We may imagine that the interaction between the anyons al- 
ters the topological liquid in which the anyons are excitations, 
producing an annulus of a different liquid. Some of the local 
fields correspond to the tunneling of a r-anyon between the in- 
ner and outer edge of the annulus. Such a process is actually 
forbidden as it would change the topological charge y. Thus, 
only fields in the trivial topological sector y = 1 are allowed 
as perturbations. This excludes ip = El^r and ij) = s' ls' r- 

Outiook Extensions to chains of anyons in topological liq- 
uids described by SU{2)k Chern-Simons theory with fc > 3 
corresponding (up to phases) to the non-abelian statistics of 
higher members of the Read-Rezayi series have been men- 
tioned below Eq. (|7]i, but their topological stability is an open 
issue. In analogy to quantum spin chains, additional inter- 
actions such as dimerization or coupling of two Fibonacci 
chains in a ladder geometry should lead to gapped quantum 
liquids. Disordered anyonic chains are currently being inves- 
tigated [I25I]. For 2D anyonic structures gapless phases of non- 
Fermi liquid type might potentially also emerge. 
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